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The theory of doped excitonic insulators is reinvestigated in light of recent experiments on hexa- 
borides. For the appropriate valley-degenerate X^^X'-g, band structure, "intra- valley" condensation 
is energetically favored. Ferromagnetism occurs upon doping due to the quenching of kinetic energy 
at the otherwise direct Erst-order excitonic insulator-metal transition. The phase diagram includes 
states of spatially inhomogeneous density and magnetization at low temperatures. 
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PACS numbers; 71.10.Ca, 71.35.-y, 75.10.Lp 

Recent dramatic measurements by Fisk et. ai. m on 
very lightly-doped divalent hexaborides (Laa^Cai-a^Bg, 
La^^Sri-j^Bg, etc.), have revived interest in the physics 
of excitonic insulators, which are coherent condensates 
particle-hole excitations that may occur in semi-metals 
with slightly overlapping or gapped conduction and va- 
lence bands. The theory of such excitonic insulators was 
developed in the 1960's and 70's (see ReL [^ for a re- 
view and original references). Band structure calcula- 
tions predict that CaBg and SrBg are indeed at the bor- 
der between very small gap semiconductors and slightly 
overlapping electron-hole metals (with gap/overlap of 
-OMeV < Eg < GMeV) |||. These materials are thus 
prime candidates for an excitonic instability. The experi- 
ments observed ferromagnetism below Tc ~ &QQK with a 
small moment maximized around x — 0.005. Early the- 
oretical work M predicted high-temperature ferromag- 
netism in doped El's, and the same connection was re- 
cently made for the hexaborides. [Q 
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FIG. 1. Schematic phase diagrams in temperature versus 
chemical potential /x and doping x supposing the EI is fully 
gapped for a; = 0. The thick dashed lines indicate first-order 
boundaries. If for x = the system is not fully gapped, a 
region of triplet EI persists at T = for small x. The order 
of the phase boundaries near the point A is sensitive to the 
detailed parameters of the model. 

In this letter the theory of ETs is reinvestigated, tak- 
ing into account the three-pocket band structure appro- 
priate to the hexaborides. A number of significant errors 
in the previous work on a two-band model are also cor- 
rected. In particular, we find that a first-order transition 
as chemical potential is varied and a consequent jump in 



the equilibrium density is crucial in the physics of ferro- 
magnetism (see Fig. 1). 

To proceed, we make the crucial physical assumption 
that the characteristic energy scales of the EI are much 
smaller than the bandwidth. A natural parameter char- 
acterizing this smallness is the ratio of the band gap Eq 
{Eg < for overlapping bands) to the bandwidth W. 
For |r| = |i?G|/T4^ ^ 1, the physics is controlled by 
small momentum-exchange processes. In this limit, the 
long-range part of the Coulomb potential dominates, and 
much of the behavior becomes insensitive to the details of 
the band structure. To make this explicit, we expand the 
electron field operators c„,cj, (a =t, i is a spin index): 



cl(r)- 



'/'.Q.+kW^L(k)- 



(1) 



Here (j)ik(r) is a Bloch function at quasi- momentum 
k with i = 1,2 for valence and conduction states, 
respectively. The integration {J^ = J ,^^K ) is re- 
stricted to spheres of radius A around the regions 
of band proximity, centered at the three wavevectors 
Qa = (tt, 0,0), (0,7r,0), (0, 0,7r) (we measure all lengths 
in units of the lattice spacing) . The "flavored" fermion 
fields (Fourier transformed back to real space) obey 
{V'iaaW^V'jb/sl'"')} = S^jSabSafiSir-r'). The effective 
mass approximation is valid for F ^ 1, since then we 
can take A <C 1. Thus H = Hq + He, with the non- 
interacting Hamiltonian 
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In the Bloch basis, H, 
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describing the long-range part of the Coulomb interac- 
tion, and 
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In the r ^ 1 limit, one finds that G can be approximated 
as momentum-independent and 0{e'^ /Q^), where Q ~ 
0(1) is a reciprocal lattice vector. Thus H' <^ He — 
0(e^/g^), and can be regarded as a small perturbation. 
The approximate "continuum" Hamiltonian He = 
Hq + He has much higher symmetry than the complete 
one. Most obviously, He conserves separately the charge 
and spin of the electrons of each flavor, which consti- 
tutes a [C/(2)] invariance. Note that He has a much 
larger U{12) symmetry, which is however strongly bro- 
ken by the inequivalent dispersion relations in Hq. As 
emphasized by Halperin and Rice g, one may visual- 
ize the system at this level as a collection of (flavored) 
positrons and electrons, which naturally have a strong 
tendency to bind into "atoms" - excitons. It is thus nat- 
ural to formally deflne an "s-wave" excitonic-insulating 
state by the presence of an off-diagonal expectation value 
(V'ioa('')V'2b/3(i")) 7^ foi' some a, 6, a,/?. Clearly, for this 
order parameter to represent a true symmetry-breaking, 
there must be no band mixing (at the points Qa) in the 
kinetic energy. Mixing is prevented for the hexaborides as 
the conduction and valence states at X belong to distinct 
{X^,X^) representations of the group of the wavevector. 
To proceed, we define local pair fields: 
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Here and in the rest of the paper we introduce Pauli 
matrices r, a acting in the band and spin spaces, respec- 
tively, and 3x3 matrices T^°' {a — 1 ... ua) acting in the 
"flavor" space. We also suppress indices wherever pos- 
sible. The excitonic order parameters are off-diagonal 
in the band space, and so utilize only the off-diagonal 
Pauli matrices with z/ = 1, 2. The appropriate T's (which 
specify irreducible representations of the cubic group - 
r;,r((,r'5,r^ for a = 1,2,3,4) are U(3) generators. In 
particular, T^^ = I/V^ and the remainder are expressed 
in terms of the standard Gell-Mann basis {Ai-'-Ag}: 
r2i = A3/V2, r22 ^ Ag/A T31 = Ai/A T32 = 

A4/V2, r33 = A4/y2, T41 = X2/V2, T42 = As/A 

T43 = \t/V2. These satisfy Tr T-^^T^^ = S'^^S"''^. Due 
to the separate U(2) invariances for each ia in ipiaa, the 
order parameters are unified for He into only two inde- 
pendent multiplets. The first contains the diagonal (in 
flavor) order parameters with ^^,^^4, with A ~ 1,2; 
the second contains the remaining off-diagonal order pa- 
rameters with A — 3,4:. 

The energetics separation of the above two multiplets is 
entirely due to the anisotropy in dispersion at the differ- 
ent Qa points. Several points suggest that this anisotropy 



generally favors the diagonal order parameters. A some- 
what physical argument is that electrons and holes that 
can best move together tend to pair more effectively. This 
suggests that states with similar group velocity at equal 
momenta (necessary for a zero-momentum condensate) 
and hence similar dispersion will tend to preferentially 
bind. Note that anistropies in the interactions also exist, 
but are much weaker than those in the dispersion. 

These heuristic considerations are supported by two 
concrete calculations. Firstly, consider approach the EI 
state by reducing Eq > in the band insulator. In 
this limit, the nature of the condensate is determined by 
the lowest-energy bound state (which reaches the chem- 
ical potential flrst). It is then necessary to compare 
the binding energies of electron-hole pairs taken from 
the same and different Kq. As a caricature of the true 
band structure, we consider a toy model with ml j^ = 
rriu ^ = JTiii^ (c.f. Eq. |3|). For an electron- hole pair 
drawn from the same X point, the effective masses match, 
and the reduced mass in the corresponding anisotropic 
hydrogen atom problem is nir = {mu/2,m±/2,m±/2). 
If the electron and hole are taken from different X 
points, one obtains instead rrir = {mj_/2,rh,rh), with 
m = TO||TOj_/(to|| + m±). Various methods (e.g. varia- 
tional) can be used to convince oneself that the "diag- 
onal" electron-hole pair is more strongly bound when- 
ever mil 7^ rn±. This is particularly clear in the limit 
TO|| ^ m±, in which case the diagonal atom has one very 
heavy reduced mass component, while all the reduced 
masses become equal to m±/2 in the off-diagonal case. 

A second argument is provided by diagrammatic tech- 
niques which apply deep into the metallic limit when 
Eg < and Coulomb effects are weak, e'^^/rn/Ea <C 1. 
It is then possible to integrate out modes and further 
reduce the momentum-space cut-off to a set of spherical 
shells of width 2 A' <C kp around the Fermi surfaces [pj . 
With the reduced cut-off, fermi-surface renormalization 
group (RG) arguments pO| or more traditional diagram- 
matics may be used to show that the only (marginally) 
relevant two-particle interactions (i.e. which lead to log- 
arithmic instabilities) are valence-conduction "ladder" 
terms in which a conduction electron is scattered to a 
valence state at a nearby point on the Fermi surface and 
vice versa. The most general allowed terms of this type 
can be written 
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where the superscript indicates the spherical shell cut- 
off, and V — \,2 sums only over the inter-band r matri- 
ces. Furthermore, to simplify the presentation we have 
kept only the most relevant s-wave component of each 
interaction channel; including the full wave-vector de- 
pendence is straightforward and leads to no signiflcant 
modiflcations. Note that with the sign convention above 



y > is an effectively "attractive" interaction. Due to 
its high symmetry, the continuum interaction Hq gives 
identical contributions to all the coupling constants, so 

yuAp/a = vb + SV^''^P/'', with Vo 
and 



27re2 ln(fcF/A')/fc|. 
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^^ — ^ilkj ^ji^k^ba ^ dc 
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Note that the corrections to the singlet interactions, 
SV^P, have a repulsive "direct" contribution absent in 
the triplet (a) channel, so that the correction terms 
favor triplet pairing, in agreement with the Hartree- 
Fock theory of Ref. ||]. Amongst the favored flavor- 
diagonal representations, there are four possible triplet 
order parameters, with A = 1,2, /i = 1,2. To de- 
cide between them requires an explicit calculation of the 
Coulomb matrix elements, and hence knowledge of the 
Bloch states. The simplest wavefunctions consistent with 
the XsjX^ symmetries are 0i = v8cos7ra;sin27r2/sin27rz, 
(f>2 — V^ sin 7ra;(cos 27rj/— cos 2ttz) for the valence and con- 
duction states at (tt, 0,0), respectively. For these wave- 
functions, one finds the most favorable correction occurs 
to the A = 2,v = 1 couphng constant, suggesting an 
EI in the Fg representation. The presence of r^ pair- 
ing indicates broken time-reversal symmetry consistent 
with a spatially-varying spin density within the unit cell. 
This result should not be taken as definitive, and might 
be modified by e.g. including higher harmonics in the 
Bloch functions. It is crucial to remember, however, that 
the splittings due to the 5V terms are only weak pertur- 
bations to the continuum contribution from He- Thus 
although a particular EI state is energetically preferred, 
the other nearly degenerate states can play important 
roles. 

Consider next a mean- field (MF) approach [^,Q. Fol- 
lowing the above reasoning, we include only the diagonal 
interactions, V'^'^p ^V"^" , with :/ = 1, 2, A = 1, 2. To 
allow analytic progress, we also assume the toy model 
dispersion above and coUinear spin polarization of all 
non-zero order parameters along the a^ axis. Define 
gap functions A^„ = f (V^Lr^'Vaa) = {Ks + "A^H/S- 
The mean-field hamiltonian neglecting 5V corrections is 

Hmf = /k Ti-MFOi-) with, 

Hmf = E (V'L Kr' - A^ - A^rl i^a. + ^ [ALl'j • 

(11) 

Because H-mf is a sum of independent terms for dif- 
ferent flavors and spin orientations, each such species 
can be studied separately. To proceed, we note that 
the hamiltonian for a single species can be mapped to 
the problem of a BCS superconductor in a Zeeman field. 
Under the particle/hole transformation ijj2aa ~^ ''t'laa: 



A]^^ + iA^^ becomes the complex BCS order parameter 
and /i the Zeeman field, if valence and conduction states 
are re-interpreted as up and down spins. This MF the- 
ory (MFT) was solved by Larkin and Ovchinnikov H and 
Fulde and Ferrell M. For small /i at low temperatures, a 
BCS state occurs, with |A| = Ao = 2tcxp{-l/NoVo) at 
T = 0, where t = vp^' is an energy cut-off and A'o is the 
density of states at the Fermi level. In the simplest MFT 
assuming a uniform order parameter, at zero tempera- 
ture a first order transition occurs at /i = /ic = Ao/\/2 
directly into the normal state. A more detailed examina- 
tion shows that an intermediate state with a non-uniform 
order parameter [|6|j^] exists in a narrow region around 
fi — fie, with a first order transition to the BCS state at 
fjL = l-i'e < fie and a possibly second order transition to 
the normal state at /i = /i" > fie- Note that the "non- 
trivial" solution used in Refs. [Q,g| is actually an unstable 
and unphysical free energy maxiuium. 

For our purposes, it is sufficient at this stage to con- 
sider only the uniform MFT - in any case the non- 
uniform solution is invalidated by the correction terms 
soon to be added to Eq. |l^. The above analysis indicates 
that for fj, = fie, there are two degenerate minima of the 
single species free energy density 
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Here (3 = l/{kBT) and the quasiparticle energies are 
^± ~ —fi±yje1. + A^A^. Since the total free energy den- 
sity is the sum /tot = T.aa /«"' ^^^^ implies a very large 
(2^^-fold) degeneracy for the full electron-hole system. 
In particular, the amplitude for each flavor and spin ori- 
entation may be chosen independently to equal or Aq. 
While two equal free energy minima occur at any first 
order transition, this large degeneracy is non-generic. 

Additional interactions can split the degeneracy in fa- 
vor of a ferromagnetic state. The search for symmetry- 
breaking terms is simplified by the fact that, because 
the different states have macroscopically different par- 
ticle occupations, only operators diagonal in flavor and 
spin have non- vanishing matrix elements in the degener- 
ate subspace. A little examination shows that the domi- 
nant perturbations are 
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A W^P term, which represents the long-wavelength 
Coulomb potential, was neglected, as appropriate for uni- 
form states provided the background charge is taken into 
account (but see below). Terms involving ^V^i/' were 



omitted from Eq. |3|, as they are important only for 
band-gap renormalization. Likewise, the 6V corrections 
directly couple the Aaa order parameters, but are negli- 
gibly small: SV < W^p/" w Vq. 

Treating the terms in Eq. O as perturbations, one 
finds that to the leading order approximation in which 
W^P = W^" = W^", the lowest energy states for the 
cubic problem comprise two sets: 8 with ferromagnetic 
(but but arbitrarily oriented) polarizations of all flavors 
(e.g. Aa+ = Ao, Aa_ = 0), and 12 with one X-point un- 
gapped, one X-point polarized, and one fully gapped (e.g. 
Ai± = 0, A2+ = Ao, Aa- = 0, Agi = Aq). The next or- 
der corrections {SW^p < SW'^'^ < 0, SW^" > 0) favor the 
aligned and fully polarized state (note that transverse 
magnetization M A A^^ 7^ is never favored) . 

To understand the behavior with chemical potential in 
more detail, we expand the MET and introduce order pa- 
rameters decoupling the W^"' and W'^'^ spin interactions 
(the W^P interaction can also be decoupled, but does not 
influence the mean-field solution in the physical parame- 
ter regime). The extended MET is obtained by replacing 
fi —>■ ^ + aha for flavor, spin a, a and adding the terms 
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^/i = I^ihl + hl + hl)-^{hih2 + h2h3 + hih3) (14) 



to the original hamiltonian density Ti.MF- Here J — 
{2TTe'^/K'^) ■> 6 ^ e^lQ'^ > 0. It is straightforward to 
solve the extended MET for S -^ J a.t zero temperature. 
One finds a ferromagnetic ground state in the range fici = 
AicVl - 2NoJ < M < Mc2 = McV(l-2AroJ)/(l-4iVoJ), 
with first order transitions at fici, Mc2 (this assumes fic2 < 
Ao). At finite temperature the MET gives the phase di- 
agram shown in Eig. 1. Note that the above behavior 
is stable to further perturbations such as the SV correc- 
tions, except insofar as to stabilize a particular (e.g. E3) 
triplet state for fi < ^^i- 

A salient feature of the MET is the first-order nature 
of the low-temperature boundary of the ferromagnetic 
phase. Taking into account long-range Coulomb inter- 
actions, macroscopic charge neutrality requires that as 
a function of doping x the system must therefore form 
an inhomogeneous (i.e. domain or labyrinthine) state 
for Q < X < Xc = 6Nofici, the minimum value in the 
ferromagnetic phase. It is straightforward to estimate 
the typical width i of ferromagnetic domains (their de- 
tailed morphology is much more difficult to determine) 
which is controlled by a competition between electro- 
static energy and surface tension a ^ NqVpAq, which 
set the profile of the electrochemical potential. One 

finds e - [a/{e^x'^{l-x/xc)^)]^^^. Note that if i be- 
comes of order the "coherence" length ^ = vf/Aq, this 
semiclassical analysis is invalid and a microscopic quan- 
tum calculation is required. A consequence of the do- 
main state is that the magnetization density is linear 
in the doping in this region, i.e. m/x is constant (and 



stays non-zero as x ^ O"*"). Similar behavior holds for 
x'^ = 6NoHc2 < X < x'^ = 12Nofj,c2, where a non- uniform 
mixture of ferromagnetic and normal domains obtains. 

Beyond the simplest MET, it is possible for the sys- 
tem to remain in an unpolarized and uniform excitonic 
state for very small doping at zero temperature. This 
can occur only if the undoped excitonic state is itin- 
erant (and therefore gapless). To study the feasibility 
of this scenario, we have considered an unnested MET 
(i.e. with different effective masses for electrons and 
holes). In this case, an intermediate-coupling excitonic 
metal phase occurs near fi — 0, with a non- vanishing 
excitonic order parameter but an incompletely gapped 
Eermi surface. A first order T = transition to an ex- 
citonic ferromagnet remains if the incommensurability is 
not too large. It is likely that a continuous excitonic 
paramagnet-ferromagnet transition is prohibited by the 
singular quasiparticle-mediated interactions, but this is- 
sue warrants further investigation. 

Many open questions remain. An enormous multitude 
of collective (pseudo-Goldstone modes of [t/(2)]'^) modes 
should exist in the \T\ <C 1 limit, mostly with small gaps. 
There are also a large number of metastable free-energy 
minima which could lead to substantial hysteresis. Non- 
trivial topological excitations are also possible. Einally, 
the role of disorder and possible trapping by dopant ions 
should be considered. The present analysis provides a 
framework for systematic studies of these issues, and 
of the applicability of the excitonic model to the hex- 
aborides. 



ACKNOWLEDGMENTS 



We thank Z. Eisk, H. R. Ott, and L. Gorkov for dis- 
cussions. 



[1] 
[2] 



D. P. Young et. ah, Nature 397, 412 (1999). 

B. I. Halperin and T. M. Rice, in Solid State Physics, 21, 

F. Seitz, D. Turnbull, and H. Ehrenreich, eds.. Academic 

Press, New York, 1968. 

S. Massidda et. al, Z. Pliys. B 102, 83 (1997); A. 

Hasegawa and A. Yanase, J. Phys. C 12, 5431 (1979). 



[7] 
[8] 
[9] 

[10] 



M. E. Zhitom irsky, T. M. Rice and V. I. Anisimov, cond- 
mat/9904330| unpublished. 

B. A. Volkov et. aL, Sov. Phys.-JETP 41, 952 (1976); 
ibid., 43, 589 (1976). 

A. I. Larkin and Yu. N. Ovcliinnikov, Sov. Phys.-JETP 
20, 762 (1965). 

P. Fulde and R. A. FerreU, Phys. Rev. 135, A550 (1964). 
L. Balents, in preparation. 

A proper treatment of the long-range Coulomb interac- 
tion in the RG approach will be deferred to Ref. |8|. 
R. Shankar, Rev. Mod. Phys. 66, 129 (1994). 



